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^ ' Abstract 

Prolate spheroidal wave functions (PSWFs) play an important role 
in various areas, from physics (e.g. wave phenomena, fluid dynamics) to 
engineering (e.g. signal processing, filter design). Even though the sig- 
nificance of PSWFs was realized at least half a century ago, and they 
I frequently occur in applications, their analytical properties have not been 

investigated as much as those of many other special functions. In particu- 
lar, despite some recent progress, the gap between asymptotic expansions 
I and numerical experience, on the one hand, and rigorously proven explicit 

bounds and estimates, on the other hand, is still rather wide. 

This paper attempts to improve the current situation. We analyze 
the differential operator associated with PSWFs, to derive fairly tight 
estimates on its eigenvalues. By combining these inequalities with a num- 
ber of standard techniques, we also obtain several other properties of the 
^ . PSFWs. The results are illustrated via numerical experiments. 

• Keywords: bandlimitcd functions, prolate spheroidal wave functions, Priifer 

I transformation 
'sf ■ Math subject classification: 33E10, 34L15, 35S30, 42C10 

■ 1 Introduction 

The principal purpose of this paper is to provide proofs for several inequalities 
•r-n ^ involving bandlimitcd functions (see Section [3] below). While some of these 

/\ . inequalities are known from "numerical experience" (see, for example, [S], |10| . 

I [TT], dH), their proofs appear to be absent in the literature. 

■ ■ ■ A function / : K — >■ R is bandlimitcd of band limit c > 0, if there exists a 

function a G L"^ [—1,1] such that 

fix) = J'^ a(t)e'-* dt. (1) 

In other words, the Fourier transform of a bandlimitcd function is compactly 
supported. While ^ defines / for all real a;, one is often interested in bandlim- 
itcd functions, whose argument is confined to an interval, e.g. — 1 < x < 1. 



*This author's research was supported in part by the AFOSR grant #FA9550-09-l-0241 
tYale University, 51 Prospect st. New Haven, CT 06511. Email: andrei.osipoviOlyale.edu. 
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Such functions are encountered in physics (wave phenomena, fluid dynamics), 
engineering (signal processing), etc. (see e.g. [2], [TH], HO])- 

About 50 years ago it was observed that the eigenfunctions of the integral 
operator Fc : L'^ [—1,1] — > [— 1, 1], defined via the formula 

Fc M (x) = J'^ V'(t)e'^-* dt, (2) 

provide a natural tool for dealing with bandlimited functions, defined on the 
interval [—1,1]. Moreover, it was observed (see, for example, [5], [T^, [12]) 
that the eigenfunctions of Fc are precisely the prolate spheroidal wave functions 
(PSWFs), well known from the mathematical physics (see, for example, [16], 
[TS^). The PSWFs are the eigenfunctions of the differential operator Lc, defined 
via the formula 

icM(x) = -^((l-x^).g(:.))+cV. (3) 

In other words, the integral operator Fc commutes with the differential opera- 
tor Lc (see [9], [H]). This property, being remarkable by itself, also plays an 
important role in both the analysis of PSWFs and the associated numerical 
algorithms (see, for example, [3], [J). 

It is perhaps surprising, however, that the analytical properties of PSWFs 
have not been investigated as thoroughly as those of several other classes of 
special functions. In particular, when one reads through the classical works 
about the PSWFs (see, for example, [5], [TU], [H], [H], [H]), one is amazed by 
the number of properties stated without rigorous proofs. Some other properties 
are only supported by analysis of an asymptotic nature; see, for example, [6], 
[7], [15], [17]. To some extent, this problem has been addressed in a number of 
recently published papers, for example, [5], [?], [5]. Still, the gap between nu- 
merical experience and asymptotic expansions, on the one hand, and rigorously 
proven explicit bounds and estimates, on the other hand, is rather wide; this 
paper offers a partial remedy for this deficiency. 

This paper is mostly devoted to the analysis of the differential operator Lc, 
defined via ([3]). In particular, several explicit bounds for the eigenvalues of 
Lc are derived. These bounds turn out to be fairly tight, and the resulting 
inequalities lead to rigorous proofs of several other properties of PSWFs. The 
analysis is illustrated through several numerical experiments. 

The analysis of the eigenvalues of the integral operator Fc, defined via 
requires tools different from those used in this paper; it will be published at 
a later date. The implications of the analysis of both Lc and Fc to numerical 
algorithms involving PSWFs are being currently investigated. 

This paper is organized as follows. In Section [21 we summarize a number of 
well known mathematical facts to be used in the rest of this paper. In Section O 
we provide a summary of the principal results of this paper. In Section |4l we 
introduce the necessary analytical apparatus and carry out the analysis. In 
Section [51 we illustrate the analysis via several numerical examples. 
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2 Mathematical and Numerical Preliminaries 



In this section, we introduce notation and summarize several facts to be used 
in the rest of the paper. 



2.1 Prolate Spheroidal Wave Functions 

In this subsection, we summarize several facts about the PSWFs. Unless stated 
otherwise, all these facts can be found in [4], [5], [7], [9], [101 . 

Given a real number c > 0, we define the operator Fc : [—1,1] — > [—1, 1] 
via the formula 

M {x) = y ^(t)e*^^* dt. (4) 

Obviously, Fc is compact. We denote its eigenvalues by Aq, Ai, . . . , A„, . . . and 
assume that they are ordered such that |A„| > |A„+i| for all natural n > 0. 
We denote by ^pn the eigenfunction corresponding to A„. In other words, the 
following identity holds for all integer n > and all real — 1 < a; < 1: 

Xni'nix) = I ^Mt)e''''' dt. (5) 

We adopt the conventiorQ that ||?An|lL2[-i.i] = 1- The following theorem de- 
scribes the eigenvalues and eigenfunctions of Fc (see [3], [5], [S]). 

Theorem 1. Suppose that c > is a real number, and that the operator Fc is 
defined via ^ above. Then, the eigenfunctions '00 j V'l > ■ • ■ of Fc are purely real, 
are orthonormal and are complete in 1,1]. The even-numbered functions 
are even, the odd-numbered ones are odd. Each function ipn has exactly n simple 
roots in (—1,1). All eigenvalues A„ of Fc are non-zero and simple; the even- 
numbered ones are purely real and the odd-numbered ones are purely imaginary; 
in particular, A„ = | A„ | . 

We define the self-adjoint operator Qc : [— 1. 1] [— 1, 1] via the for- 

mula 

^ r 1 / X ^ sin (c (x — t)) , , , 

Qc M (x) = - / , fit) dt. (6) 

TT X - t 

Clearly, if we denote by J' : L^(]R) L'^{M.) the unitary Fourier transform, then 

Qc M (x) = xi-ui] (x) ■ [xi-c.c] iO ■ ^ if] iO] (^), (7) 

where X[-a,a] : K — >■ R is the characteristic function of the interval [—a, a], 
defined via the formula 

I otherwise. 



This convention agrees with that of [1], [5] and differs from that of [9]. 
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for all real x. In other words, Qc represents low-passing followed by time- 
limiting. Qc relates to Fc, defined via (|H), by 

Qc^^-f:- f„ (9) 

and the eigenvalues /i„ of Qn satisfy the identity 

Mn = ^-|A«|', (10) 

for all integer n > 0. Moreover, Qc has the same eigenfunctions ■0n as F^- In 
other words, 

, , , 1 /"^ sin(c(a; — i)) , , , , ,,,, 

t^ni^nix) = - / ' Mt) dt, (11) 

J-1 X-t 

for all integer n > and all — 1 < x < 1. Also, Qc is closely related to the 
operator Pc : L^(R) — > L^(M), defined via the formula 

„ r w N 1 /""^ sin (cix — t)) , . , 
Pe M (x) = - / ^ V ¥'(t) dt, (12) 

which is a widely known orthogonal projection onto the space of functions of 
band limit c > on the real line R. 

The following theorem about the eigenvalues of the operator Qc, defined 
via dl]), can be traced back to [7]: 

Theorem 2. Suppose that c > and < a < 1 are positive real numbers, and 
that the operator Qc '■ 1, 1] — > L'^ [—1, 1] is defined via © above. Suppose 
also that the integer N{c, a) is the number of the eigenvalues Hn of Qc that are 
greater than a. In other words, 

iV(c, a) = max {fc = 1,2,... : ^i,._i>0}. (13) 

Then, 

N{c,a) = -c+ ( ^log^— ^ ) logc + 0(logc). (14) 



According to (HH), there are about 2c/7r eigenvalues whose absolute value is 
close to one, order of log c eigenvalues that decay exponentially, and the rest of 
them are very close to zero. 

The eigenfunctions ^„ of Qc turn out to be the PSWFs, well known from 
classical mathematical physics (see, for example, [12], [IS])- The following the- 
orem, proved in a more general form in [12j . formalizes this statement. 

Theorem 3. For any c > 0, there exists a strictly increasing unbounded se- 
quence of positive numbers Xo < Xi < • ■ • such that, for each integer n > 0, the 
differential equation 

(1 - x^) ij/'ix) - 2x ■ i;'{x) + {xn - c^x^) i^{x) = (15) 

has a solution that is continuous on [—1,1]. Moreover, all such solutions are 
constant multiples of the eigenfunction ipn of Fc, defined via ^ above. 
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The following theorem provides lower and upper bounds on %„ of Theorem [3] 
(see, for example, [1], [5], [TO]). 

Theorem 4. For all real c > and all natural n> 0, 

n{n + l) <Xn <n{n + l) + c^. (16) 

The following result provides an upper bound on "0^(1) (see [5])- 
Theorem 5. For all c > and all natural n> 0, 

rnil)<n+l. (17) 



2.2 Elliptic Integrals 

In this subsection, we summarize several facts about elliptic integrals. These 
facts can be found, for example, in section 8.1 in '8 , and in [21^. 

The incomplete elliptic integrals of the first and second kind are defined, 
respectively, by the formulae 



(18) 



F{y,k)= / , ^ , 
Jo vl-fc sin i 

E{y, k) = r Vl- sin^ t dt, (19) 

^0 

where < 2/ < 7''/2 and < k < 1. By performing the substitution x = sint, 
we can write (fTSl) and (ITOl) as 

r'"*^^^ dx 

^^^'^^~io ^{l-x^){l-k^x^)' ^^^^ 

Eiy,k) = J^ \l^—^dx. (21) 

The complete elliptic integrals of the first and second kind are defined, respec- 
tively, by the formulae 

Fik)^F(l,k)= r (22) 

/TT N , 

E{k) = E \^-,kj = J Vl-k^sin^tdt, (23) 
for all < fc < 1. Moreover, 

E (v/T^) = 1 + (^-1 + log(2) - . fc2 + O . log(fc)) . (24) 



5 



2.3 Oscillation Properties of Second Order ODEs 

In this subsection, we state several well known facts from the general theory of 
second order ordinary differential equations (see e.g. p^). 

The following two theorems appear in Section 3.6 of [I] in a slightly different 
form. 

Theorem 6 (distance between roots). Suppose that h{t) is a solution of the 
ODE 

y"{t) + Q{t)y{t)=Q. (25) 
Suppose also that x < y are two consecutive roots of hit), and that 

< Q{t) < B^, (26) 

for all X < t < y . Then, 

Theorem 7. Suppose that a < b are real numbers, and that (7 : (a, &) — > R is a 
continuous monotone function. Suppose also that y{t) is a solution of the ODE 

y"{t)+g{t)-y{t)^0, (28) 

in the interval (a, 6). Suppose furthermore that 

ti<t2<t3<... (29) 

are consecutive roots ofy(t). If g is non- decreasing, then 

t2-t,>t3~t2>t4~t3> .... (30) 

// g is non-increasing, then 

t2-ti<t3-t2<ti-t3< .... (31) 



The following theorem is a special case of Theorem 6.2 from Section 3.6 in [T]: 

Theorem 8. Suppose that gi, g2 are continuous functions, and that, for all real 
t in the interval (a, 6), the inequality gi{t) < g2{t) holds. Suppose also that 
01, 02 : (a, fe) — > R are real-valued functions, and that 

(l>';it)-^-g,it)-Mt)^0, 

02W+.92(i)-02(<)=O, (32) 
for all a < t < b. Then, 02 has a root between every two consecutive roots of 

01. 

Corollary 1. Suppose that the functions 0i,02 are those of Theorem\^ above. 
Suppose also that 

0i(io)-02(to), 0'i(to) = 02(io), (33) 
for some a < to < b. Then, 02 has at least as many roots in (tg, b) as 0i. 
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Proof. Due to Theorem [51 we only need to show that if ti is the minimal root of 
in (to, b), then there exists a root of 02 in (^Oj ^i)- By contradiction, suppose 
that this is not the case. In addition, without loss of generality, suppose that 
4'i[t) 1 4'2[t) are positive in (to,ti). Then, due to (15^ . 

0102 - 02^1 = (52 - 51) '?^i'/'2, (34) 

and hence 

< / (32(5) -gi(s))0l(s)02(s)rfs 

J to 

= [0'i(s)02(s)-0i(s)0^(s)]*: 

= 0'i(ti)02(ti) <0, (35) 
which is a contradiction. ■ 



2.4 Priifer Transformations 

In this subsection, we describe the classical Priifer transformation of a second 
order ODE (see e.g. [I],[12]). Also, we describe a modification of Priifer trans- 
formation, introduced in [3] and used in the rest of the paper. 
Suppose that in the second order ODE 

^^[p{t)u\t))+q{t)u(t)=Q (36) 

the variable t varies over some interval / in which p and q are continuously 
differentiable and have no roots. We define the function 6* : / — > M via 

P{t)u'{t) 

-— = 7 t tan6'(t , 37 

u{t) 

where 7 : / — ^ M is an arbitrary positive continuously differentiable function. 
The function 0{t) satisfies, for all t in /, 

*).,:2„,,> l(t) /'T'(<)^ ™(2«(0) 



One can observe that if u'{t) = for t E I, then by p7)) 

9[t) — kir, k is integer. (39) 

Similarly, if u{i) = Q for i E I, then 

0{i) = {k + 1/2) TT, k is integer. (40) 

The choice j{t) = 1 in ([57|) gives rise to the classical Priifer transformation (see 
e.g. section 4.2 in [T]). 

In [3], the choice 7(t) — q{t)p{t) is suggested and shown to be more con- 
venient numerically in several applications. In this paper, this choice also leads 
to a more convenient analytical tool than the classical Priifer transformation. 

Writing in the form of ([55)1 yields 

p{t)^l-t\ q{t)^Xn-cH\ (41) 
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for \t\ < min {^/Xn/c, l}. The equation ([37]) admits the form 



P^MM = V^(i)^tan^?W, (42) 

Ipnit) 



which imphes that 



where m{t) is an integer determined for aU t by an arbitrary choice at some 
t = to (the role of ■nm{t) in is to enforce the continuity of at the roots of 
7/i„). The first order ODE ^ admits the form (see [5]. [^1 

0'(<) = -/(i) + sin(20(i))<t), (44) 

where the functions /, v are defined, respectively, via the formulae 



and 



1 p(0g^(0 + g(0y(t) 1 / i ^ 



3 Summary 

In this section, we summarize some of the properties of prolate spheroidal wave 
functions (PSWFs), proved in Section H] The PSWFs and the related notation 
were introduced in Section l2.ll Throughout this section, the band limit c > 
is assumed to be a positive real number. 

Many properties of the PSWF %pn depend on whether the eigenvalue x„ of 
the ODE fTS]) is greater than or less than c^. The following simple relation 
between c, n and Xn is proven in Theorem [TS] in Section 14.1.21 

Proposition 1. Suppose that n > 2 is a non-negative integer. 

* If n < (2c/7r) — 1, then Xn < ■ 

* If n'> {2c/ t:), then Xn > (? ■ 

* // (2c/7r) — 1 < n < (2c/7r), then either inequality is possible. 



In the following proposition, we describe the location of "special points" 
(roots of V'n, roots of tp'^^, turning points of the ODE (fT5)) l that depends on 
whether Xn > or Xn < ■ It is proven in Lemma [1] in Section 14.1.11 and is 
illustrated in Figures [U [5] 

Proposition 2. Suppose that n > 2 is a positive integer. Suppose also that 
ti < ■ ■ ■ < tn are the roots of ipn in (—1, 1), cind xi < ■ ■ ■ < Xn-i are the roots 
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of ip'„ in {ti,tn). Suppose furthermore that the real number Xn is defined via the 
formula 



j maximal root of ijj^-^ in (—1, 1), if Xn < c^, 
ll, ifxn>c^- 



(47) 



Then, 



'Xn yXn /<o\ 
< -Xn <ti<Xi<t2<--- < tn-1 < tn < Xn < • (48) 



C 

^2 



In particular, if Xn < c , then 

tn <Xn< ^ < 1, (49) 
C 

and ip'n has n + 1 roots in the interval (—1, 1); and, if Xn > c^, then 

tn<Xn^l< (50) 
C 

and ip'n has n — 1 roots in the interval (—1, 1). 

The following two inequalities improve the inequality (|16|) in Section 12.11 
Their proof can be found in Theorems I9I10I in Section 14.1.21 This is one of the 
principal analytical results of this paper. The inequalities ([?T|) . ([5^ below are 
illustrated in Tables [H H H H 

Proposition 3. Suppose that n > 2 is a positive integer. Suppose also that tn 
and T are the maximal roots of ipn and ipn interval (—1, 1), respectively. 

If Xn > c^, then 



TT 



2 /•*" I xn - cH^ , ^2 I xn - cH^ 



// Xn < (? , then 



TT 



1 + 3/ ,r-^dt<n<- d^f^dt. (52) 



A'^oie that (|5ip ami (|52p differ only in the range of integration on their right-hand 
sides. 

In the following proposition, we simplify the inequality (jSip in Proposition[3l 
It is proven in Theorem [T51 and Corollary [3] in Section B. 1.31 

Proposition 4. Suppose that n > 2 is a positive integer, and that Xn > c^. 
Then, 



TT 







^ /x;:-i;f^ ) <n + 3, (53) 



where the function E : [0, 1] — !■ R is defined via (j23p m Section 
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The following proposition is an immediate consequence of Proposition It 
is proven in Theorem [^D] in Section 14.1.31 and is illustrated in Figures [31 131 

Proposition 5. Suppose that n > 2 is a positive integer such that n > Icj-n , 
and that the function f : [0,oo) — M zs defined via the formula 

f{x) = -l+ Vt + cos2(6') de. (54) 

Suppose also that the function H : [0,oo) — )■ R is the inverse of f , in other 
words, 

y = f{H{y)) = -1 + / v/^(2/)+cos2(0) dO, (55) 
Jo 

for all real y > 0. Then, 




In the following proposition, we describe a relation between Xn a-nd the 
maximal root i„ of V'n in (~ljl)j by providing lower and upper bounds on 
1 — t„ in terms of Xn and c. It is proven in Theorem [T71 [TH] in Section 14.1.31 

Proposition 6. Suppose that n > 2 is a positive integer, and that Xn > (? ■ 
Suppose also that i„ is the maximal root of ipn in the interval (—1, 1). Then, 

^2/8 



X„ - C2 + V(Xn - C2)2 + (^c/2) 



2 



4^2 

< , (57) 

Xn - c2 + ^{Xn - c2)2 + (47rc)2 



The following proposition is a special case of Proposition [51 It is proven in 
Theorem [23l in Section 14.1.31 and is illustrated in Figure [5l 

Proposition 7. Suppose that c > IO/tt. Suppose also that n > 2 is a positive 
integer, and that 

n > - + 1 + i • log(c). (58) 

TT 4 

Suppose furthermore that tn is the maximal root of ipn in the interval (—1,1). 
Then, 

1 . 27r2 

<l~tn< J. (59) 



8 ■ (1 + V2) X« - c2 Xn - c- 



In the following proposition, we provide yet another upper bound on Xn in 
terms of n. Its proof can be found in Theorem [T^l in Section I4.1.3[ and it is 
illustrated in Tables [SI [51 and Figure [31 



10 



Proposition 8. Suppose that n > 2 is a positive integer, and that Xn > <? ■ 
Then, 

Xn<(|(n+1))'. (60) 



We observe that, for sufBciently large n, the inequality (j60p is even weaker 
than ([TB)) . On the other hand, (|60l) can be useful for n near 2c/7r, as iUustrated 
in Tables El El 

The following proposition summarizes Theorem I 111 in Section r4. 1.21 and The- 
orems [H [HI in Section SH It is illustrated in Tables |S1 EHHIHl 

Proposition 9. Suppose that n > 2 is a positive integer, and that Xn > <? ■ 
Suppose also that —l<ti<t2<---<tn<l are the roots oftpn in the interval 
(— 1, 1). Suppose furthermore that the functions f,v are defined, respectively, via 
(|l5l),(IMl) 171 Section\K^ Then: 

• For each integer (n + l)/2 < i < n — 1, i.e. for each integer i such that 
< < tn, 

f{t.+,)+v{t,+,)/2<'''^'-''^jk)- ^^^^ 

• For each integer (n+\)/2<i<n— \, i.e. for each integer i such that 

< < tn, 

ti+1 — ti > ti+2 — ti+1 > ■ ■ ■ > tn ~ ^n-l- (62) 

• For all integer j = 1, . . . ,n — 1, 

The following proposition summarizes Theorem 1161 in Section [4.1.3l 

Proposition 10. Suppose that n > 2 is an integer, and that Xn < — c\f2. 
Suppose also that — 1 < ti < i2 < ■ ■ • < in < 1 are the roots of ^„ in the interval 
(—1,1). Then, 

ti+1 — ti < ti+2 — ti+1 < ■ ■ ■ < tn — tn-1, (64) 

for all integer (n+l)/2 < i < n— 1, i.e. for all integer i such that < ti < tn. 

The following proposition summarizes Theorem [5l in Section 12.11 and Theo- 
rem |2l in Section 1121 



Proposition 11. Suppose that n > is a non-negative integer, and that Xn > 
(? . Then, 

i<^^(l)<n+i. (65) 
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The following proposition is illustrated in Figures [TJ [21 It is proven in The- 
orem [55] in Section 14.1.31 

Proposition 12. Suppose that n > is a non-negative integer, and that x,y 
are two arbitrary extremum points of ipn in (—1,1). If \x\ < \y\, then 

\Mx)\ < \My)\ ■ (66) 

//, in addition, Xn > , then 

\Mx)\<\My)\<\M)\- (67) 



4 Analytical Apparatus 

The purpose of this section is to provide the analytical apparatus to be used in 
the rest of the paper, as well as to prove the results summarized in Section [31 

4.1 Oscillation Properties of PS WFs 

In this subsection, we prove several facts about the distance between consecutive 
roots of PSWFs (O and find a more subtle relationship between n and Xn ifTSj) 
than the one given by (fTH)) . Throughout this subsection c > is a positive real 
number and n is a non-negative integer. The principal results of this subsection 
are Theorems [SI [TUl [H and [T3 

4.1.1 Special Points of i/)„ 

We refer to the roots of ■i/'n, the roots of ■0^ and the turning points of the ODE 
([T5|) as " special points" . Some of them play an important role in the subsequent 
analysis. These points are introduced in the following definition. 

Definition 1 (Special points). Suppose that n > 2 is a positive integer. We 
define 

• ti < t2 < ■ ■ ■ < tn to be the roots of tpn in (—1, 1), 

• xi < ■ ■ ■ < a;„_i to be the roots of ip'„ in (ii,i„), 

• Xn via the formula 

I maximal root of ip'^ in (—1, 1), 



This definition will be used throughout all of Section [H The relative location 
of some of the special points depends on whether Xn > <? or Xn < <? ■ This is 
illustrated in Figures [H [21 and is described by the following lemma. 



(68) 
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Lemma 1 (Special points). Suppose that n > 2 is a positive integer. Suppose 
also that ti < ■ ■ ■ < tn and xi < ■ ■ ■ < Xn are those of Definition]^ Then, 



< -Xn < ti < Xi < t2 < ■ ■ ■ < tn-1 < tn < X„ < . (69) 



C 

^2 



In particular, if Xn < c , then 

tn <Xn< ^ < 1, (70) 
C 

and ip'^ has n + 1 roots in the interval (—1, 1); and, if Xn > , then 

tn <Xn = \< (71) 
C 

and ip'n has n — 1 roots in the interval (—1, 1). 
Proof. Without loss of generality, we assume that 

^n(l) > 0. (72) 

Obviously, ((72)) implies that 

iP'nitn) > 0. (73) 

Suppose first that Xn < c^- Then, due to the ODE (|T5]) in Section [2Jl 

V';(l) = ^^^-V'n(l)<0. (74) 

We combine ([15]) and (|73p to obtain 

<(*») = — ^•V';(U>0. (75) 

In addition, we combine ([S5|) . ([75)) . (fM)) to conclude that the maximal root x„ 
of V'^i in (—1, 1) satisfies 

tn<Xn< 1. (76) 
Moreover, ([7^ implies that, for any root x of ip'^ in (tn, 1), 



C(^) = - ^" 'y -^„(x)<0. (77) 

1 — 



We combine ([72]), ([761), dZZj) with ([15]) to obtain 

1 ^n '0n(*^n) 



< 0, (78) 



which implies both (l69l) and ([70)) . In addition, we combine ([T5)). ([70| and ([77)) 
to conclude that x„ is the only root of ip'n between tn and 1. Thus, ip'^ indeed 
has n + 1 roots in (—1,1). 
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Suppose now that Xn > c^. We combine (ITS]) and ([7^ to obtain 

^;(l) = ^^^-^n(l)>0. (79) 
If i„ < a; < 1 is a root of ■0^^, then 

C {^)=- ^" ~ • V-n (X) < 0, (80) 

therefore i/'^ can have at most one root in (t„, 1). We combine this observation 
with (|68l) . ([75]) . ([75]) and (IM)) to conchide that, in fact, t/;^ has no roots in 
{tn, 1), and hence both and (|7T]) hold. In particular, -0^ has n — 1 roots in 
(-1,1). ■ 

4.1.2 A Sharper Inequality for Xn 

In this subsection, we use the modified Priifer transformation (see Section [Z4)) 
to analyze the relationship between n, c and Xn- In particular, this analysis 
yields fairly tight lower and upper bounds on Xn in terms of c and n. These 
bounds are described in Theorems I9I10I below. These theorems are not only 
one of the principal results of this paper, but are also subsequently used in the 
proofs of Theorems HH [H [H [HI [H 

We start with developing the required analytical machinery. In the following 
lemma, we describe several properties of the modified Priifer transformation (see 
Section [231) , applied to the prolate differential equation 

Lemma 2. Suppose that n > 2 is a positive integer. Suppose also that the 
numbers ti, ... ,tn and Xi, . . . ,Xn are those of Definition{l\in Section \4-.l.l\ and 
that the function 6 : [— a:„,x„] — >■ IR. zs defined via the formula 



6{t) 



(i — • TT, if t = ti for some 1 < « < n, 

atan ■ + m(t) • n, z/V«(t)y^O, 

(81) 

where m{t) is the number of the roots of tpn in the interval (—l,t). Then, 9 has 
the following properties: 

• 9 is continuously differentiable in the interval [— a;„,a;„]. 

• 9 satisfies, for all — a;„ < t < Xn, the differential equation 

9'{t) = f{t) + v{t)-8m{29{t)), (82) 



where the functions f,v are defined, respectively, via (j45p . (j46p in Sec- 
tion\2^. 



• for each integer < j < 2n, there is a unique solution to the equation 

^W-J-f, (83) 
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for the unknown t in [—Xn,Xn]- More specifically, 

(84) 

^, (85) 
(86) 

for each i — 1, . . . ^n. In particular, 9(xn) — n ■ i:. 

Proof. We combine in Lemma [T] with (ISTI) to conclude that 9 is well defined 
for all —Xn < t < Xn, where Xn is given via (|55| in Definition [T] Obviously, 
is continuous, and the identities ((84)) . ((85)) . ((86)) follow immediately from the 
combination of Lemma [T] and ((8T|) . In addition, satisfies the ODE (j82p in 
{-Xn,Xn) due to dM]), dH]) , (gll) in Section!^ 

Finally, to establish the uniqueness of the solution to the equation ([55]) . we 
make the following observation. Due to (j81L for any point t in (— a;„,a;„), the 
value 9{t) is an integer multiple of tt/2 if and only if t is either a root of V'n or 
a root of i/'^ . We conclude the proof by combining this observation with , 
(IHl and jHll). ■ 

Remark 1. We observe that, due to ((84|) . ([85)) . ([M]) . /or all i = 1, . . . ,n, 

sin(26l(t,)) = sin(26'(a;,)) = 0, (87) 

where ii, . . . , Xi, . . . , x„ are t/iose o/ Definition]^ in Section \4-.l.l\ and 6 is 
that of Lemma [H This observation will play an important role in the analysis 
of the ODE (|82p throughout the rest of this section. 

In the following lemma, we prove that 9 of Lemma [5] is monotonically in- 
creasing. 

Lemma 3. Suppose that n > 2 is a positive integer. Suppose also that the real 
number Xn and the function 9 : [— a;„,a;„] — !■ M are those of Lemma\^ above. 
Then, 9 is strictly increasing in [— a;„,a;„], in other words, 

9'{t) > 0, (88) 

for all —Xn <t < Xn- 

Proof. We first prove that 

s(7)(">°- W 

for —Xn < t < Xn, where the functions f,v are defined, respectively, via (1451) . 



9{-Xn) = 0, 
9(Xi) = i • TT, 
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in Section We differentiate vj J with respect to t to obtain 

'p'q + q'p fp\' fp'q + q'p\' 



Apq 



j3/2pl/2 



q ^p ^ 



{p"q + 2p'q' + q"p) q^'^p^'^ 



g-5/2p-3/2 



\q'p + \p'q ] {p'q + q'p) - pq {p"q + 2p'q' + q"p) 



\p^q'f + \q\p'f^q'pp"-p\q" 



>0, 



(90) 



since, due to (|4T|) . 

p{t) > 0, p"{t) = -2 < 0, q{t) > 0, = -2c2 < 0. (91) 

We now proceed to prove ([55)1 for < f < x„. Suppose that, by contradiction, 
there exists < x < Xn such that 



9'ix) < 0. 

Combined with (|82p in Lemma [2] above, (j92p impHes that 

, ^ ^'(a^) . ^^ /(a;) + i;(a;) ■ sin(2g(x)) 
1 + • sm(2t^(x)) = j-—^ < 0, 



(92) 



(93) 



and, in particular, that 

sin(26'(a;)) < 0. (94) 
Due to Lemma [2] above, there exists an integer (n + l)/2 < i < n such that 

i - • vr < 9{x) < i ■ TT. (95) 

Moreover, due to (gS]), dM]), (gl]), (|M1), there exists a point y such 

that 

0<ii<a;<j/<Xi< a:„, (96) 

and also 

e{x) = e{y), e'{y)>o. (97) 

for otherwise ([551) would be impossible. We combine ([5^ and ([W)) to obtain 



1 + ^ . sin(2.(x)) = /(^) + -(^);^^in(2%)) ^ ^ ^ 



/(y) 



f{y) 



fiv) 



(98) 



in contradiction to ([89|), ([93]) and ([Ml)- This concludes the proof of ([88|) for 
Q <t < Xn- For —Xn < t < 0, the identity ([55)) follows now from the symmetry 
considerations. ■ 
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The right-hand side of the ODE of Lemma contains a monotone term 
and an oscillatory term. In the following lemma, we study the integrals of the 
oscillatory term between various special points, introduced in Definition [1] in 
Section HXn 

Lemma 4. Suppose that n > 2 is an integer. Suppose also that the real numbers 
ti < ■ ■ ■ < tn and xi < ■ ■ ■ < Xn are those of Definition]^ in Section \4.1.1\ and 
the function 6 : [— a;„,x„] — >■ M is that of Lemma\^ above. Suppose furthermore 
that the function v is defined via (j46|) in Section \2.4\ Then, 



v{t) ■ sin(26l(t)) dt > 0, (99) 

v{t) ■ sin(26l(t)) dt < 0, (100) 

ti+i 

v{t) ■ sin(26l(t)) dt < 0, (101) 

for all integer (n — l)/2 < i < n — 1, i.e. for all integer i such that < Xi < Xn- 
Note that the integral in (jlOip is the sum of the integrals in (j99p and (jlOOp . 

Proof. Suppose that z is a positive integer such that (?i — l)/2 < i < n — 1. 
Suppose also that the function s : [0, ti • tt] — > [— a;„,x„] is the inverse of 9. In 
other words, for all < rj < n ■ tt, 

0{s{v)) - V- (102) 

Using dMl), (ESI), §B in Lemma d we expand the left-hand side of ^ to 
obtain 



v{t) ■ sin(26l(t)) dt = 

e{u+i) 

v{s{ri)) ■ sin(277) • s'iji) drj ■ 

(^.) 

y[s{r^)) . sin(2r7) drj 



f{s{r^))+v{s{r,))-sin{2r,) 
u(s(r/ -1- i ■ tt)) • sin(2?7) dry 
f{s{v + i ■ tt)) -f v{s{r] -I- i • tt)) • sin(277) ' 



(103) 



from which ([M)) readily follows due to in Section and ([55| in Lemma [31 
By the same token, we expand the left-hand side of (jlOOD to obtain 

v{t) ■ s\n{2e{t)) dt = 

ti+i 

'^'■'^^^"^ v(s{r])) ■ sin(277) dr] 
(i+i/2y^ fisiv)) + Hsiv)) ■ sin(2r7) 

v{s{v + (» + 1/2) ■ tt)) ■ sin(27y) drj 

/o f{s{v + {i + 1/2) ■7r))-v{s{7j + {i + 1/2) ■7r))-sm{2r])' ^ ' 
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which, combined with in Section and ([BS]) in Lemma [31 imphes (jlOOp . 
Finally, for all < 77 < 7r/2, 

sin(2r7) ^ sin(277) ^^^^^ 



+ + 1/2) • tt)) - sin(27?) (//w)(s(7/ + z • vr)) + sin(2r7) ' 

since the function f /v is decreasing due to (15^)) in the proof of Lemma |31 The 
inequality (jlOip now follows from the combination of (I103D , p04D and (|105p . ■ 

We are now ready to prove one of the principal results of this paper. It is 
illustrated in Tables [D [H HIH 

Theorem 9. Suppose that n>2 is a positive integer. If Xn > <? , then 



n<\l\\^^dt. (106) 



TT 







Ifx-n. < c^, then 



TT 



2 r IXn-cH^ 



where T is the maximal root of in (—1,1). Note that p06p and (|107p differ 
only in the range of integration on their right-hand sides. 

Proof. Suppose that the real numbers 

-1 < -Xn < ti < Xi < t2 < ■ ■ ■ < tn-1 < Xn^l < tn < Xn < 1 (108) 

are those of Definition [1] in Section [4.1.11 and the function 9 : [— x„,a;„] — ?► M 
is that of Lemma [5] above. Suppose also that the functions /, v are defined, 
respectively, via (HSl) , in Section 12.41 If n is even, then we combine (j82l) , 
((55)) . (I5S)) in Lemma [2] with (fTUTj) in Lemma H to obtain 



TT = / e\t) dt^ f{t) dt+ ^ v{t) ■ sm{20{t)) dt 

< / " fit) dt. (109) 



If n is odd, then we combine (|H2]), (gS]), dSl]) in Lemma H with ([TOIIl) . pUT]) in 
Lemma m to obtain 



n 

2'" 



9'{t)dt^ / f{t)dt + 

(Ti + l)/2 -^0 



/ v{t) ■ sin(20(t)) + 2J / ■ sin(26'(t)) 

< / " /(i) dt. (110) 
Jo 

We combine (I109P and (IllOp with (HSl) in Lemma [T] to conclude both (|106p and 
(|TU71) . ■ 
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To prove Theorem [TUl we need to develop a number of technical tools. In 
the following two lemmas, we describe several properties of the equation f(t) = 
v{t) in the unknown t, where f,v are defined, respectively, via (|^5|) . pS)) in 
Section!^ 

Lemma 5. Suppose that n > is a non-negative integer. Suppose also that the 
functions /, v are defined, respectively, via (|45l) , (j46p in Section \2.4\ Suppose 
furthermore that the real number Xn is that of Definition [7] in Section \4.1.1\ 
Then, there exists a unique point t in the interval (0, a;„) such that 

fii)^vii). (Ill) 



Proof We observe that, due to (gS]) , (gS]) in Section 

vit) 



for all < i < x„ . Moreover, 



4^ = 0, hm 44=-- (113) 
/(O) t^.„.t<.„ f(t) ^ ' 

We combine ((89l) in the proof of Lemma [3] with (jll2p and (|113p to conclude 
both existence and uniqueness of the solution to the equation f{t) = v{t) in the 
unknown < t < Xn- B 

Lemma 6. Suppose that n > 2 is a positive integer. Suppose also that the real 
number Xn and the function 9 : [~XmXn] — > M are those of Lemma\^ above. 
Suppose furthermore that the point < t < Xn is that of Lemma\^ above. Then, 



TT < e{i) <n-TT. (114) 



Proof. Suppose that the point < x < x„ is defined via the formula 

"-J)-^), (115) 

where denotes the inverse of 9. By contradiction, suppose that (|114p does 
not hold. In other words, 

0<i <x. (116) 

It follows from the combination of Lemma [5l (j89p in the proof of Lemma [Sj and 
(jliep . that f{x) < v{x). On the other hand, due to (|S2|) in Lemma [5] and (jllSp . 

9'{x) = f{x)+v{x) ■ sin(26l(a;)) 

= f{x) + v{x) ■ sin (2mT - |) = /(^) - v{x) < 0, (117) 

in contradiction to ([55)) in Lemma [31 ■ 
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In the following three lemmas, we study some of the properties of the ratio 
//?;, where f,v are defined, respectively, via (^5]) . pS)) in Section [ 



Lemma 7. Suppose that n > is a non-negative integer, and that the functions 
f,v are defined, respectively, via pS)) . PS|) in Section Then, for all real 

< t < 1, 



dt\v 



(t) = ht{a) ■ fit), (118) 



where the real number a > is defined via the formula 

Xn 

a = — , 

and, for all < t < I, the function ht : (0, oo) — > K is defined via the formula 

At^ + (2a - 6) • + (4 - 8a) • + 2a • (a + 1) 
<2.(l + a-2i2) 

Moreover, for all real < t < minjy^, 1}, 



(119) 



-ti, ^ i^z,u, - • -r v-i - ou,; • t -r • V" T J-/ . 



4(i)w>^/W^ (121) 



Proof. The identity (jll8p is obtained from (|^5|). via straightforward alge- 
braic manipulations. To establish (I12ip . it suffices to show that, for a fixed 
< t < 1, 

3 

inf |/it(a) -.t^ <a<oo]>-. (122) 

a 2 

We start with observing that, for all < i < 1, 

2 

lim /if (a) = 6, lim ht{a) — —. (123) 

a^t^, a>t^ a->-oo 

Then, we differentiate ht{a), given via (|120p . with respect to a to obtain 

It follows from (jl23p . (I124p . that if < at < oo is a local extremum of /if (a), 
then 

-6t4 + 9f2_i 
= > * ' (125) 

which is possible if and only if 1 > > 1/7. Then we substitute at, given via 
(fT25)) . into ([T20| to obtain 

— + 14/2 _ 1 

Mt,at)^ \ (126) 

It is trivial to verify that 

inf |/i(t, at) : ^ < t < l| = Jim^^ /i(/, at) = | (127) 
Now follows from the combination of p^ . p^ . and P?7)) . ■ 
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Lemma 8. Suppose that n > 2 is a positive integer, and that t„ is the maximal 
zero ofTpn is the interval (—1, 1). Suppose also that the real number Zq is defined 
via the formula 

Zo = ~ 0.4591. (128) 

Then, for all < t < tn, 

v{t) < fit) ■ Zo, (129) 
where the functions f,v are defined, respectively, via (|45|) . (|46p in Section \2.4\ 

Proof. Due to ([59)1 in the proof of Lemma [31 the fmiction f /v decreases mono- 
tonically in the interval (0,i„), and therefore, to prove (|129p . it suffices to show 
that 

Suppose that the point i is that of Lemma[ni Suppose also that the real number 
Xn and the function 9 : [—Xn,Xn\ — M are those of Lemma [51 Suppose further- 
more that the function s : [O.u-tt] — !■ [— x„,a;„] is the inverse of 9. In other 
words, for all < ?/ < n • tt, 

9{s{rj)) = r?. (131) 
We combine Lemma [H Lemma [21 Lemma [5l Lemma [51 and Lemma [3 to obtain 



* I = ^ (~i) (^(^)) 

^dt\ v)^^' fKv)) + v{s{v)) ■ sin(27y) > 



(n-l/4)7r 



dt 



(-1) {s{v}) dv 



(«-i/2)7r fisin)) + vis{7j)) ■ sin(2?7) 



> 



(-V^)- i (-{) i^iv)) dy n 3^3n 
/(„-i/2). fisiv)) ^ 4 ■ 2 - 8 ' ^ ^ 

which implies (|130p . ■ 

Lemma 9. Suppose that n > 2 and (n + I) /2 < i < n — 1 are positive integers. 
Suppose also that the real number Xn and the function 6 : [— x„,x„] — ^ M are 
those of Lemma\^ Suppose furthermore that Q < 5 < 7r/4 is a real number, and 
that the real number Zs is defined via the formula 

-1 



Zs 



(133) 



2 V 4 1 + Zq- sin(2(5) ^ 

where Zq is defined via p28p in Lemma \E\ above. Then, 

vit) < fit) ■ Zs, (134) 

for all < t < s ((i + 1/2) ■ n — S), where the functions /, v are defined, respec- 
tively, via (|45p . (|46p in Section \2.4\ and the function s : [0, n • tt] — )■ [— x„,a;„] 
is the inverse of 9. 
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Proof. Suppose that the point tg is defined via the formula 

ts = s {{i + 1/2) -TT-d). (135) 

Due to (j89D in the proof of Lemma [21 the function f/v decreases monotonically 
in the interval (0,^^), and therefore to prove (|134p it suffices to show that 



v{ts) Zs 2 V4 1 + • sin(2(5) , 

We observe that, due to Lemma |3l 

< sin(26'(t)) < sin(25), (137) 

for all ts <t < s((i + l/2)7r). We combine ([HSt . with Lemmad LemmaH 

Lemma [51 Lemma [7] and Lemma [S] to obtain 

fits) /(g((z + l/2)^)) ^ 
v{ts) w(s((i + l/2)7r)) 

d f f\ _ /-(^-i/a)- i (-() (s(77)) dry 

V ^) '^^ " Ji^-l/2).-s fisiv)) + y(Mv))-M2v) ^ 



(-1/2). ^(-{)(s(r?)) 



dri 

> 



[^-l,2).-S fisiv)) 1 + ivlf)isirj)) ■ sin(25) 

3 1 

2"^-l + Zo.sin(2<5)- ^^^^^ 

We combine (jl38l) with (I129p in Lemma[5]to obtain p36L which, in turn, implies 
(fml) . ■ 

In the following two lemmas, we estimate the rate of decay of the ratio f/v 
and its relationship with 9 of the ODE ([5^ in Lemma [5J 

Lemma 10. Suppose that n > 2 and {n+ 1)/2 < i < n — 1 are positive integers. 
Suppose also that the real number Xn and the function 9 : [— a;„,x„] — )■ ]R are 
those of Lemma\^ Suppose furthermore that < d < 7r/4 is a real number. 
Then, 

- ) (s(i7r - S)) - (-] (s(i7r -6 + n/2)) > 2 ■ sm{26), (139) 



where the functions f,v are defined, respectively, via (|3S]), (H^ in Section 
and the function s : [0, n • tt] — > [— a;„, Xn] is the inverse of 9. 

Proof. We observe that, due to Lemma [5] and Lemma [21 

sin(26'(t)) > 0, (140) 

for all s(i7r) < t < s{iTT — S + 71/2). We combine (jl40l) with Lemma[51 Lemma[21 
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Lemma [HI Lemma [7] and Lemma 1^] to obtain 

fisjiTT)) fis{i7r-S + 7T/2)) _ 
w(s(i7r)) w(s(«7r — 5 + 7r/2)) 

"-"^^^iHli!^ >^.f^-.).^, (141) 

f{s{v)) l + {v/f){s{r^j} 2 \2 J 1 + Zs' ^ ^ 

where Zg is defined via p33p in LemmalHl We also observe that, due to Lemma[2] 
and Lemma [3l 

sin(26l(i)) < 0, (142) 

for all s(i7r — (5) < t < s{iTr). We combine (|142l) with Lemma HJ Lemma [31 
Lemma [5] and Lemma [7] to obtain 

fisjin-S)) fisjin)) ^ 
v{s{i'K — 5)) v{s{iT:)) 

d r f \ _ n Tt (~{) (sW) d-n ^ 

s(t^-S) dt\ v) Ji^_s /(s(?7)) + v{s{t])) ■ sin(2r7) 

Next, suppose that the function ft, : [0, 7r/4] — > R is defined via the formula 

M^) = ^-(f-^)-^ + ^-^-2.sin(2^), (144) 

where Zs is defined via (|133p in Lemma [SI One can easily verify that 

mm{h{5) : 0<5<n/4:}> (145) 
S 2d 

and, in particular, that h(6) > for all < S < n/A. We combine ((TiT|) . 
(flii)) and (fTI5| to obtain ([1351) . ■ 

Lemma 11. Suppose that n > 2 and (n + l)/2<i<n—l are positive integers. 
Suppose also that the real number Xn and the function 6 : [— — > IR are 
those of Lemma\^ Suppose furthermore that Q < 5 < 7r/4 is a real number. 
Then, 

^ (s(i7r + 5- Ti/2)) - {s{in + 6j) > 2 ■ sin(2J), (146) 

where the functions /, v are defined, respectively, via (1451) , (146^ in Section 
and the function s : [0, n • tt] — > [— a;„, is the inverse of 9. 
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Proof. We observe that, due to Lemma [5] and Lemma [31 

sm{29(t)) > 0, (147) 

for all s(i7r) < t < s{iTi + 5). We combine (I147P with Lemma [H Lemma |31 
Lemma [SI Lemma [7| and Lemma [Hi to obtain 

w(s(i7r)) u(s(i7r + (5)) 

s{^.+s) ^ / ^ i (-|) (g(^)) dy ^ 

V J in /(s(?7)) + w(s(?7)) • sin(277) 

i (-1) i^M) 3 ^ 

/(s(77)) ■ 1 + (V/)(s(ry)) 2'l + Z,' ^ ^ 

where Zg is defined via (jl33p in Lemma[9l We also observe that, due to Lemma[2] 
and Lemma [3l 

sin(26'(t)) < 0, (149) 

for all s{iTT + S — n/2) <t< s{in). We combine (|149p with Lemma [U Lemma [31 
Lemma [5] and Lemma [3 to obtain 

f{s{in + 5-Ti/2)) fisjin)) _ 
visjin + S — tt/2)) v{s(i7r)) 

^(-) d_ / f\ _ ft (-{) i^i^)) ^ 

f{s{r])) +v(s{r])) ■ sin(2?7) 
ft (-1) (^(^)) dy 3 /vr X 

»7r+5-7r/2 f[s{v)) 2 \2 J 

Obviously, for all < S < n/i, 

3 6 3 /TT x 3 /TT \ 1 3 , 

>o- o-'5 •— -7-+o-'5- (151) 



2 1 + Za- 2 V2 / 2 V2 / 1 + 2 

We combine (fT48|) . (fTSOj) . ([TsT]) with (fT44)) . ([T45)) in the proof of Lemma [10] to 
obtain (ITi5|) . ■ 

In the following lemma, we analyze the integral of the oscillatory part of 
the right-hand side of the ODE between consecutive roots of ipn- This 
lemma can be viewed as an extention of Lemma [4l and is used in the proof of 
Theorem [TUl below. 

Lemma 12. Suppose that n > 2 is an integer, —l < ti < t2 < ■ ■ ■ < tn < 1 are 
the roots of ipn in the interval (—1, 1), and xi < ■ ■ ■ < a;„_i are the roots of 
in the interval {ti,tn)- Suppose also, that the real number Xn and the function 
6 : [—Xn,Xn] — 7> K are those of Lemma\^ above. Suppose furthermore that the 
function v is defined via (j46p in Section \2.4\ Then, 

ti+i 

v{t) ■ sin(26l(t)) dt > 0, (152) 

ti 

for all integer {n + l)/2 < i < n—1, i.e. for all integer i such that < ti < tn. 
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Proof. Suppose that i is a positive integer such that {n ~ l)/2 < i < n — 1. 
Suppose also that the function s : [0,n • tt] — > a;„,a;„] is the inverse of 9. In 
other words, for all < 77 < n • tt, 

9{s{v)) = V- (153) 
Due to p04p in the proof of Lemma 2] above, 



v{t) ■ sin{20{t)) dt = 

v{s{iTT + f] - tt/2)) ■ sm{2r]) dri 
Q /(s(i7r + r) — 7r/2)) — v{s{iiT + rj — 7r/2)) • sin(277) 



(154) 



We proceed to compare the integrand in (jl54p to the integrand in (|103p in the 
proof of LemmaHl First, for all < 77 < 7r/4, 

< ^, / ^^ , ■ ^ . (155) 



(//w)(s(i7r + 77 - 7r/2)) - sin(277) (//w)(s(i7r + 77)) + sin(277) 

due to p46p in Lemma [TT] Moreover, for all 7r/4 < 77 < 7r/2, we substitute 
5 — Ti /2 ~ rj to obtain 



< 



{f/v){s{iTT + 77 - tt/2)) - sin(2r7) {flv){s[m - (5)) - sin(2(5) 

{f/v){s{m -5 + ^/2)) +sixi{25) " (//v)(s(i^ + ?7))+sin(2r7)' ^^^'^^ 

due to (I139P in Lemma [TUl We combine (|103l) in the proof of Lemma |3] with 
(fT5i)l . (fT55)l . (|T5B1) to obtain (fT^ . ■ 

The following theorem is a counterpart of Theorem |9] above. It is illustrated 
in Tables mil El m 

Theorem 10. Suppose that n>2 is a positive integer. Suppose also that i„ is 
the maximal root of ipn in (—1, 1)- Then, 



TT 



2 r Xn-cH^ 



1 + z/ \l \^t^ dt<n. (157) 



Proof. Suppose that the real numbers 

-1 < "Xn < ti < Xi < t2 < ■ ■ ■ < tn-1 < Xn-1 < t„ < X„ < 1 (158) 

and the function : [— a;„,a;„] — > M are those of Lemma [2] above. Suppose also 
that the functions /, w are defined, respectively, via (|^5|) . p5)) in Section [531 If 
n is odd, then we combine ([5^ . ([55|) . in Lemma [5] with (jl52p in Lemma [T^ to 
obtain 

— :^ — • TT = 



/" 9'{t)dt^ f{t)dt+ f v{t) ■ sm{29{t)) dt 

> f{t)dt. (159) 
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If n is even, then we combine ([5^ . (|55|) . ([BS]) in Lemma with (|99p in Lemma [3] 
and (|152p in Lemma [T^ to obtain 



n ^ 1 

Z ■ TT 



e'{t) dt^ f{t) dt + 



/ v{t) ■ sm{2e{t)) dt+ / v{t) ■ sin(26l(t)) 

> f" fit)dt. (160) 



We combine ([75^1) and (ITM| to conchide (ITFf)) . ■ 

Corollary 2. Suppose that n > 2 is a positive integer, and that Xn > <? ■ 
Suppose also that tn is the maximal root of ipn in the interval (—1, 1). Then, 

I + -VX^ ■ E ( asm{t„) , <n< -y/x^ - E [ , (161) 

where E{y, k) and E{k) are defined, respectively, via (|19p and (j23p in Sec- 
tion[ 



Proof. It follows immediately from (I?!]) , (jl06l) in Theorem [S] and (|157p in The- 
orem [ini ■ 

The following theorem, illustrated in Tables [71 [51 provides upper and lower 
bounds on the distance between consecutive roots of ^„ inside (—1, 1). 

Theorem 11. Suppose that n > 2 is a positive integer, and that Xn > ■ 
Suppose also that — 1 < ti < <2 < ■ ■ ■ < < 1 the roots o/^„ in the interval 
(— 1, 1). Suppose furthermore that the functions f,v are defined, respectively, via 
gg),® in Section\2^\ Then, 

<t,+i-t,< (162) 



f(t,+,) + v{t.,+,)/2 ^ fit,)' 

for all integer {n + l)/2 < i < n~l, i.e. for all integer i such that < ti < i„. 

Proof. Suppose that the function : [— 1,1] ^ M is that of Lemma [21 We 
observe that / is increasing in (0,1) due to (|l5)) in Section and combine 
this observation with , (1551) , in Lemma [H and (|152l) in Lemma [T^l to obtain 

0'{t) dt = [ * f{t) dt + [ ^ v{t) ■ sm{29{t)) dt 



> ^ fit) dt > it,+i ~ t,) ■ fit,), (163) 

which implies the right-hand side of (jl62[) . As in Lemma [21 suppose that Xi is 
the zero of -0^1 in the interval (ti,ti+i). We combine (15^ . (1551 . in LemmalU 
and dM]), (jTUO)) in Lemma [1 to obtain 

/ ' fit) dt> [ ^ 9' it) dt = Tr= 9' it) dt> [ fit) dt. (164) 
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Since / is increasing in {ti,ti+i) due to (^5]) in Section the inequality (I164p 
implies that 

Xi - U > U+i - Xi. (165) 

Moreover, we observe that v is also increasing in (0, 1). We combine this ob- 
servation with ([TM|) . (EH), dSi]), in Lemma [2] and ([Ml), ([TUIIl) in Lemma H to 
obtain 

vr = ^ ^ 6''(t) < ^ ^ /(t) dt + J ^ v{t) ■ sm{29{t)) dt 

< {ti+i - U) ■ f{U+i) + {U+i - x^) ■ v{ti+i) 

< {U+i - U) ■ fiU+i) + . y(^t,+i), (166) 

which implies the left-hand side of (|162p . ■ 

The following theorem is a direct consequence of Theorem [9] above. 
Theorem 12. Suppose that n > 2 is a positive integer. If n > 2c/tt, then 

Xn > C^. 

Proof. Suppose that Xn < c^, and T is the maximal root of in (0, 1), as in 
Theorem [5] above. Then, 



2 Ixn - c2i2 2c /y„/c2 - 1^ , 2c 2c , , 



TT Jq V 1 — ^2 TT Jo V 1 — ^2 TT TT 

due to (fTUTl) in Theorem [3 ■ 

4.1.3 A Certain Transformation of the Prolate ODE 

In this subsection, we analyze the oscillation properties of ipn via transforming 
the ODE (|15p into a second-order linear ODE without the first-order term. The 
following lemma is the principal technical tool of this subsection. 

Lemma 13. Suppose that n > is a non-negative integer. Suppose also that 
that the functions 5'„, Qn ■ (^1, 1) — > K are defined, respectively, via the formu- 
lae 



= Mt) ■ Vl-t^ (168) 

and 



l-t2 



(1-^2) 



for —1 < t < 1. Then, 

Kit)+Qnit)-^nit)^0, (170) 

for all -1 <t <1. 
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Proof. We differentiate ^„ with respect to t to obtain 



Then, using (|17ip . we differentiate with respect to t to obtain 

2t 



(171) 



Kit) =^:it)VT^-Kit) 



2t 



1 -i2 



1 



Vi^ 
1 



(l-t^).C(i)-2t-CW-^ 



1 



(172) 



V 1-^' (i2-i)2y 

We observe that (fTTO)) follows from ([T7^ . ■ 

In the next theorem, we provide an upper bound on Xn in terms of n. The 
results of the corresponding numerical experiments are reported in Tables [SI [HI 

Theorem 13. Suppose that n > 2 is a positive integer, and that Xn > <? ■ Then, 

2 



n + 1 



(173) 



Proof. Suppose that the functions Qn : (—1, 1) — )■ M are those of Lemma fT3l 
above. We observe that, since Xn > , 



Quit) > Xn + 1, 



(174) 



for — 1 < t < 1. Suppose now that t„ is the maximal root of i/'n in (^Ij !)■ We 
combine p74p with p70p in Lemma \W\ above and Theorem |S1 Corollary [1] in 
Section [2.31 to obtain the inequality 



t„>l- 



Vx^TTT 

Then, we combine (jl75p with Theorem [TOl above to obtain 



(175) 



TT 







1 -t2 



TT 



> i + ^Vx;;> i + -vx;:( 1 / — 

TT TT V VXn + 1 



> -Vx^- 1' 

TT 



(176) 



which implies (jl73p . 

The following theorem is a consequence of the proof of Theorem 1131 
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Theorem 14. Suppose that n > 2 is a positive integer, and that Xn > ■ 
Suppose also that ii < • • • < t„ are the roots of ipn in (—1, !)• Then, 



for all j — 1,2, ... ,n — 1. 



Proof. The inequality (I177p follows from the combination of p74l) in the proof of 
Theorem [T3l (IT70| in Lemma [13] and Theorem [3 Corollary [1] in Section [221 ■ 

The following theorem extends Theorem [T^j in Section 14.1.21 

Theorem 15. Suppose that n > 2 is a positive integer. 

*Ifn< {2c/tt) — 1, then Xn < c^- 

* If n> {2c/ Ti), then Xn > (? ■ 

• // {2c/ t:) ~ 1 < n < {2c/ t:), then either inequality is possible. 

Proof. Suppose that Xn > c^, and that the functions ^„,(3ra : (—1, 1) — >■ K are 
those of Lemma [T51 above. Suppose also that ti < ■ ■ ■ < tn are the roots of tpn 
in (—1, 1). We observe that, due to (11691) in Lemma [T51 

Q„(i)^c^ + f^ + ^^^>c^ (178) 



We combine (|178p with (jl70p in Lemma [T3| above and Theorem [6| in Section [2731 
to conclude that 

tj+i-tj<^, (179) 

for all J = 1, . . . , n — 1, and, moreover, 

l-t„<-. (180) 
c 

We combine (jl79l) with (jl80|) to obtain the inequality 

2 (l - -) < 2tn = tn -ti< {n - 1) -, (181) 



which implies that 



n>-c-l. (182) 

TT 



We conclude the proof by combining TheoremlT^in Section l4. 1 .2l with (|182p . ■ 

The following theorem is yet another application of Lemma [T^ above. 

Theorem 16. Suppose that n > 2 is a positive integer. Suppose also that 
— 1 < ti < t2 < ■ ■ ■ < tn < 1 are the roots ofipn in the interval (—1, 1). Suppose 
furthermore that i is an integer such that Q < ti < tn, i.e. (n + l)/2<i<n— 1. 
If Xn > c^, then 

ti+l — ti > ti+2 — ti+1 > ■ ■ ■ > tn — tn-l- (183) 
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If Xn < — c\/2, then 

ti+l ~ ti < U+2 — ti+l < ■ ■ ■ < tn — tn~l- (184) 



Proof. Suppose that the functions Q„ : (— 1, 1) — > M are those of Lemma [T51 
above. If x„ > c^, then, due to (jl69p in Lemma fTSl 

is obviously a monotonicahy increasing function. We combine this observation 
with prni) of Lemma inland ((50)) of Theorem [7] in Section [Ol to conclude 
Suppose now that 

X„ < - cV2. (186) 
Suppose also that the function P„ : (1, oo) — > M is defined via the formula 



Pn{y)^Qn Ul-^j =y^ + {Xn~C^)-y + c\ (187) 



for 1 < y < CXI. Obviously, 

n (f\ = p ( 
Suppose also that yo is defined via the formula 



Qn{t)^Pn(-^\. (188) 



1 

l-(%/36r/c)^ C^-Xn 

We combine ^Ml) . (|W1) and (ITM| to conclude that, for 1 < y < yo, 

P;(y) = 2y- {c^ - Xn) < 2yo (c' - Xn) = "J"' ~ < 0. (190) 

- Xn 

Moreover, due to ([T87| . ([T89| . ([T90| . 

P„(y) > P„(2/o) = f > 0, (191) 

for all 1 < y < yo- We combine p^ . p^ . p^ . and p^ to con- 

clude that Qn is monotonicahy decreasing and strictly positive in the interval 
(0, \/Xjnlc)- We combine this observation with (|5T|) of Theorem[7]in Section 1^751 
((70)) of Lemma m and pTCT of Lemma [H] to conclude p^ . ■ 

Remark 2. Numerical experiments confirm that there exist real c > and 
integer n > such that — c^/2 < Xn < o.nd neither of (I183p . p84p is true. 



In the following theorem, we provide an upper bound on 1 — t„, where t„ is 
the maximal root of i/'n in the interval (—1, 1). 
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Theorem 17. Suppose that n > 2 is a positive integer, and that Xn > ■ 
Suppose also that i„ is the maximal root of ipn in the interval (—1, 1). Then, 

^■{l-t^f + ^^-^-{l-t,,)<ir\ (192) 

Moreover, 

47r2 

l-tn< , (193) 

Xn - C2 + V(Xn - <-^y + (4^C)2 

Proof. Suppose that the functions Qn : (—1, 1) — )• M are those of Lemma [T3l 
above. Since Xn > c^, the function Qn is monotonicahy increasing, i.e. 

Qn{tn) < Q{t), (194) 

for aWtn <t <\. We consider the solution tpn of the ODE 

^'^{t)+Qn{tn)-Vn{t)=0, (195) 

with the initial conditions 

f{tn) = -^nitn) = 0, ¥>'(t„) - K{in)- (196) 

The function has a root j/„ given via the formula 

Suppose, by contradiction, that j/n < 1- Then, due to the combination of p70p 
of Lemma [131 Theorem |8l Corollary [1] in Section 12.31 and (jl94p above, has 
a root in the interval (t„,2/„), in contradiction to (I168p . Therefore, 

a/ Qn{tn) 

We rewrite as 

(l-t„)' •Q«(tn) < (199) 
and plug p69p into (jl99p to obtain the inequality 

^^■(i-^")^ + Tr£-(^-^") + (rTu^<-^' 

which immediately yields (jl92p . Since 1 — 1„ is positive, (j200p implies that 1 — 1„ 
is bounded from above by the maximal root Xmax of the quadratic equation 

c2 . + . X - TT^ = 0, (201) 

given via the formula 

^ • ( V(X« - c2)2 + 16^V - (x„ - c^)) 



■-max — 

16^2^2 



x« - c2 + v/(Xn - c^)' + 167rV 

which implies (|193l) . 



(202) 
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The following theorem uses Theorem [T7] to simplify the inequalities (I106P in 
Theorem El and p?7|) in Theorem [TUl in Section BX^ 



Theorem 18. Suppose that n > 2 is a positive integer, and that Xn > (? ■ 
Suppose also that tn is the maximal root of ipn iri the interval (—1, 1). Then, 



2 r IXn cH^ ^ ^203) 



Moreover, 



2 f Y„.~cH 



,2+2 



n<- I \l ^ _ ^3 dt<n + 3. (204) 



Proof. We observe that, for i„ < i < 1, 
^" ^ "^^^ < ^" ~ = ^ ~ = c2(l - i„) + (205) 

l+t l + tn l+tn l+t„ ^ ' 1 + tn ^ ' 

We combine (j205p with (I192p in Theorem [T71 to obtain the inequality 



Xn - cH^ 



l+t 1-t 



< ^— , (206) 



valid for i„ < t < 1. We conclude from (|206p that 



1 IXn - cH^ ^ ^ r dt 



dt < • / —== = —== ■ 2Vl - tn = 27r, (207) 



which implies (j203l) . The inequality (I204p follows from the combination of (|203p , 
(pMI in Theorem Hand p37)) in Theorem [TU] in Section iXl ■ 

Corollary 3. Suppose that n> 2 is a positive integer, and that Xn > ■ Then, 
n<-^,- e{-^P=] <n + 3, (208) 



where E{k) is defined via (|23p in Section \2.Si 

Proof. The inequality (|208l) follows immediately from the combination of ([25)1 
in Section Holland (pOil) in Theorem flE above. ■ 

The following theorem extends Theorem [17] above by providing a lower 
bound on 1 — t„, where t„ is the maximal root of ipn in the interval (—1,1). 

Theorem 19. Suppose that n > 2 is a positive integer, and that Xn > <? ■ 
Suppose also that t„ is the maximal root of ipn in the interval (—1, 1). Then, 

, ""'^^ <l-tn. (209) 

Xn - C2 + ^{Xn - (?Y + (7rc/2)2 
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Proof. Wc combine the inequalities (|106p in TheoremlHland (|157p in Theorem [TUl 
in Section 14.1.21 to conclude that 



K?/";^'''- (-0) 



TT 



We combine (|210p with (|205l) in the proof of Theorem [TBI above to obtain 



1 < -, , , h C^(l - in) 



< -v/c2(l~i„)2 + (x„-c2).(l-i„)- (211) 



We rewrite (|21ip as 



c'(l - t„f + ixn - c") ■ (1 - t„) - ^ > 0. (212) 

16 

Since 1 — t„ is positive, (|212l) implies that 1 — i„ it is bounded from below by 
the maximal root Xmax of the quadratic equation 

2 2 

c -x -\ a; - — = 0, (213) 

given via the formula 

a^max = ^ • (v/(Xn-c2)2+ ^2^2/4 _ (^^ _ ^2)^ 



2 



^2c2 



Xn - C2 + A/(Xn-c2)2+7r2c2/4 

which implies (j209l) . 



(214) 



The following theorem is a direct consequence of Theorem[TSl It is illustrated 
in Figures O SI 

Theorem 20. Suppose that n >2 is a positive integer such that n > 2c/n, and 
that the function f : [0,oo) K is defined via the formula 

f{x) = -l+l V^ + cos2 (61) de. (215) 

Suppose also that the function H : [0, oo) — > M is the inverse of f , in other 
words, 

2/ = /(iJ(y)) = ~1 + / V^f(y)+cos2(0) de, (216) 
Jo 

for all real y > 0. Then, 

( niT \ Yn — c2 / niT 37r \ , ^ 

^(2^-^)<^<^U-l+2^ ■ (2^^) 
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Proof. Obviously, the function /, defined via (|215p . is monotonically increasing. 
Moreover, /(O) = 0, and 

lim f{x) = oo. (218) 

Therefore, H{y) is well defined for all y > 0, and, moreover, the function H 
is monotonically increasing. This observation, combined with Theorems 1151 1181 
above, implies the inequality (j217p . ■ 

In the following theorem, we provide a simple lower bound on H, defined 
via (1^ in Theorem [501 



Theorem 21. Suppose that the function H : [0, oo) — )• M is defined via (|216p 
in Theorem \2(A Then, 

s<h(^- -Xog^^y (219) 

for all real < s < 5. 

Proof. The proof of (|219p is straightforward, elementary, and is based on ([24]) 
in Section [521 it will be omitted. The correctness of Theorem has also been 
validated numerically. ■ 

Remark 3. Numerical experiments by the author indicate that the relative error 
of the lower bound in (|219p is below 0.07 for all < s < 5; moreover, this error 
grows roughly linearly with s to ~ 0.0085 for all < s < 0.1. 

In the following theorem, we provide a lower bound on Xn for certain values 
of n. 

Theorem 22. Suppose that a is a real number, and that 

< a < 5c. (220) 

Suppose also that n > 2 is a positive integer, and that 

2c ck , / 16ec\ , 

71 > - + — • log . (221) 

TT In \ a J 

Then, 

Xn>c^ + ac. (222) 



Proof. Suppose that the function H : [Q,oo) — > R is defined via (j216l) in The- 
orem [201 It was observed in the proof of Theorem [201 that H is monotonically 
increasing. We combine this observation with (j220p . (I22ip and Theorem to 
conclude that 




(223) 



Thus (15^ follows from the combination of ([2^ and Theorem [501 ■ 
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In the following theorem, we provide upper and lower bounds on 1 — t„, 
where tn is the maximal root of ipn in the interval (—1,1), in terms of Xn — ■ 
This theorem is illustrated in Figure [5l 



Theorem 23. Suppose that 



c> — . (224) 

TT 



Suppose also that n > 2 is a positive integer, and that 

n> — + l + i •log(c). (225) 

TT 4 

Suppose furthermore that t„ is the maximal root of -ipn in the interval (—1,1). 
Then, 

X« > + I • c, (226) 



and also, 



1 

<l-tn< ^. (227) 



• (1 + V2) Xn - C2 Xn- C 



Proof. We combine (HH]), (1^^ and Theorem to obtain (P^ . Then, we 
combine ([^^ with Theorems [T71 [191 to obtain (P?7)) . ■ 

4.2 Growth Properties of PSWFs 

In this subsection, we establish several bounds on l-^nl and \tpn\. Throughout 
this subsection c > is a fixed positive real number. The principal results of 
this subsection are Theorems [Ml [25l [551 The following lemma is a technical 
tool to be used in the rest of this subsection. 

Lemma 14. Suppose that n > is a non-negative integer, and that the func- 
tions p, g : M — R are defined via (|4ip in Section \2.4\ Suppose also that the 
functions Q,Q : (0, min {y^Xn/c, l}) — > M are defined, respectively, via the for- 
mulae 

Qit) = i^lit) + ^ . mt)f - m + ^^'''^'t'f^' (228) 
q{t) Xn - cH^ 

and 

Q{t) = p{t) ■ q{t) ■ Q{t) 

= (1 - t^) ■ - ch^) ■ ^Plit) + (1 - t') ■ {yj'Jt)f) . (229) 

Then, Q is increasing in the interval (O, min {y'xn/c, l}), and Q is decreasing 
in the interval (O, min {y/xn/c, l}) • 
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Proof. We differentiate Q, defined via (|228p . with respect to t to obtain 



2 



Xn - CH 

Due to (HSl) in Section EH 



^•^^ ■i^'^it)-^). (230) 



for all —1 < t < 1. We substitute (|23ip into (|230l) and carry out straightforward 
algebraic manipulations to obtain 

Q'W = _^^,2^2)2 • iXn + c'- 2ch') ■ mt)f . (232) 

Obviously, for all < i < min { ^Jxn/ c, 1 } , 

Xn + c^ - 2ch^ > 0. (233) 

We combine p32p with p33l) to conclude that 

Q'{t) > 0, (234) 

for all < t < min {y^^/c, l}. Then, we differentiate Q, defined via (|229L 
with respect to t to obtain 

Q'it) = -2t- [ixn ch^) ■ + (1 - 1') ■ iKit)f) 

+ (1 - 1^) ■ {-2ch ■ ^l{t) + 2-{xn~ ch^) ■ Mt) ■ iUt) 

-2t ■ ii^'Mf + 2-{l-t')- ^'„{t) ■ ^y:{t)) . (235) 

We substitute (|23ip into (j235p and carry out straightforward algebraic manip- 
ulations to obtain 

Q'{t) = -2t ■ [xn + - 2c2t2) . ^l{t). (236) 



We combine (|233|) with (I236P to conclude that 

Q'{t) < 0, (237) 

for all < t < min {y/x^/c,l}. We combine (I^M)) and ((257)) to finish the 
proof. ■ 

In the following theorem, we establish a lower bound on |?/'ti(1)|- 

Theorem 24 (bound on |^„(1)|). Suppose that Xn > ■ Then, 

|V„(1)| > (238) 
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Proof. Suppose that the function Q : [—1,1] — > R is defined via (|228l) in 
Lemma 1141 Then, Q is increasing in (0,1), and is continuous in [—1,1] (see 
Lemma [HI and Theorem [3] in Section [^TT|) . Therefore, 

i^l{t)<Q{t)<Q{l)=^l{l), (239) 

for all real < i < 1. Due to Theorem [1] in Section [^TTl 



1 1-1 

2 f4-\ j+ ^ / „/.2 /I \ J+ _ „/,2 / 



i'mdtK / ^^(1) (it = V'^(l), (240) 
Jo 



which implies (|238l) . 



The following theorem describes some of the properties of the extrema of ipn 
in (-1,1). 

Theorem 25. Suppose that n > is a non-negative integer, and that x,y are 
two arbitrary extremum points of %l^n in (—1,1). // < then 

\Mx)\<\My)\. (241) 

//, in addition, Xn > , then 

|^«(a;)| < |V«(y)| < |^n(l)|. (242) 



Proof. We observe that \'ipn\ is even in (—1,1), and combine this observation 
with the fact that the function Q : [— 1, 1] — > M, defined via (|228p . is increasing 
in (0, 1) due to Lemma [HI ■ 

In the following theorem, we provide an upper bound on the reciprocal of 
iV'nl (if n is even) or [■(/'^J (if n is odd) at zero. 

Theorem 26. Suppose that x„ > c^. If n is even, then 



l^n(0)| - V C- 

// n is odd, then 



^ <4.,/n.§. (243) 



1 ArT 

Proof. Since Xn > c^, the inequality 

i^l{t)<i'l{l)<n+^, (245) 
holds due to Theorem [5] in Section [2?T] and Theorem [25l above. Therefore, 

t i'l{t)dt<l + -^<^. (246) 
Ji-i/8n 8 16n 16 
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Combined with the orthonormahty of ijjn , this yields the inequaUty 



1-1/8)1 



Since 



dx 



l-l/8ri 



1 3 5 

2 " 16 ~ 16' 



{i~x^y 2 1 



1 , x + l 
4 1 — a; 



(247) 



(248) 



it follows that 



/' 

Jo 



1 



dx 



{1-x^f 



1 - l/8n 



2 1 - (1 - l/8n)2 4 



1, 2-l/8n 



l/8n 



1 8n(8n-l) 1 



2 16n - 1 
4n + n. < 5n. 



log (16ri - 1) < 



(249) 



Suppose that the functions Q{t),Q{t) are defined for —1 < t < 1, respectively, 
via the formulae (|228p . (I229p in Lemma [Til in Section 11?^ We apply Lemma [Til 
with to — and < t < 1 to obtain 



0(0) =g(o)-p(o)-g(o)-g(o) 



> Q(t) = 



(t2-l)(V;;(i))^ 



(c2-t2-;t„) 

>cV^(0(l-t^) (x„/c^-i^)>cV^(i)(l-t^) 
It follows from ([2171) . (jMSl) and ((250)) that 



(250) 



5n.g(0)-^>Q(0)-^ / 



which, in turn, implies that 



„ f'-'/''' dx 



{i-x^y 



> 



l-l/8ri 



V4(0 dt> — 



1 / 1 C 

< Ion • 



Q(o) 



16' 
(251) 



(252) 



If n is even, then ip'n{0) = 0, also, if n is odd, then '0n(O) = 0. Combined with 
(HnH), this observation yields both (|M5| and (l^Til) . ■ 



5 Numerical Results 

In this section, we illustrate the analysis of Section 2) via several numerical 
experiments. All the calculations were implemented in FORTRAN (the Lahey 
95 LINUX version) and were carried out in double precision. The algorithms 
for the evaluation of PSWFs and their eigenvalues were based on [T]. 
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We illustrate Lemma [T] in Figures [U [51 via plotting i/)„ with Xn < and 
Xn > respectively. The relations ([70)1 and ([7T|) hold for the functions in 
Figures m [21 respectively. Theorem holds in both cases, that is, the absolute 
value of local extrema of V'n(i) increases as t grows from to 1. On the other 
hand, (j242p holds only for the function plotted in Figure [21 as expected. 



v„(t) 




Figure 1: The function ipnit) for c = 20 and n — 9. Since Xn ~ 325.42 < c^, 
the location of the special points is according to (j70p of Lemma [H The points 
y/Xn/c « 0.90197 and 1 are marked with asterisks. Compare to Figure\^ 

In Tables H H [1 we illustrate Theorems [H [TO] in the case of Xn > ■ The 
band limit c > is fixed per table and chosen to be equal to 10, 100 and 1000, 
respectively. The first two columns contain n and the ratio Xn/c^- The third 
and fourth column contain the upper and lower bound on n defined, respectively, 
via (|106p in Theorem [S] and (|157p in Theorem [TOl i.e. 

Below(n) = 1 + - / \ — — dt = 1 + ■ E asin (i„) , , 

TT Jo V 1 - TT V VX^J 

Above(n) = ^ /' J ^" ~ f/ dt = ■ E (^) , (253) 

where E denote the elliptical integrals of Section 12. 2[ and t„ is the maximal 
root of f/'n iu (—1,1) (see also (1161^ ). The fifth and sixth columns contain the 
relative errors of these bounds. The first row corresponds to the minimal n 
for which Xn > ■ We observe that, for a fixed c, the bounds become more 
accurate as n grows. Also, for n = [2c/7r] + 1 the accuracy improves as c grows. 
Moreover, the lower bound is always more accurate than the upper bound. 
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v„(t) 

8i ' ■ ■ 




Figure 2: The function ipnit) for c — 20 and n = 14. Since Xn ~ 437.36 > c^, 
the location of the special points is according to (j7ip of Lemma[l\ The points 1 
o-'fT'd ypXjnjc ~ 1.0457 are marked with asterisks. Compare to Figure{^ 



n 


/ 2 
Xn/C- 


Above (n) 


Below (n) 


Ahove{n)—n 


n— Below(n) 


n 


n 


6 





10104EH 


hOl 


0.65036E+01 


0.59568E+01 


0.83927E-01 


0.71987E-02 


10 





16310EH 


hOl 


0.10498E+02 


0.99600E+01 


0.49826E-01 


0.39974E-02 


15 





29137EH 


hOl 


0.15494E+02 


0.14963E+02 


0.32940E-01 


0.24599E-02 


20 





47078E+01 


0.20495E+02 


0.19964E+Q2 


0.24737E-01 


0.17952E-02 


25 





70050EH 


hOl 


0.25496E+02 


0.24965E+02 


0.19820E-01 


0.14066E-02 


30 





98035EH 


hOl 


0.30496E+02 


0.29965E+02 


0.16538E-01 


0.11533E-02 


35 





13103EH 


h02 


0.35497E+02 


0.34966E+02 


0.14189E-01 


0.97596E-03 


40 





16902EH 


h02 


0.40497E+02 


0.39966E+02 


0.12425E-01 


0.84521E-03 


45 





21202EH 


h02 


0.45497E+02 


0.44966E+02 


0.11052E-01 


0.74500E-03 



Table 1: Illustration of Theorems\^ \1(A with c — 10. The quantities Above(n) 
and Below(n) are defined by (j253p . 
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n 




Above(ri) 


Below(n) 


Abovc(n)— n 


n — Bclow(n) 


n 


n 


64 


0.10066E+01 


0.64590E+02 


0.63964E4 


-02 


0.92169E-02 


G.56216E-03 


70 


0.10668E+01 


0.70513E+02 


0.69971E4 


-02 


0.73216E-G2 


0.40732E-03 


75 


0.11290E+01 


0.75505E+02 


0.74971E4 


-02 


0.67341E-02 


0.38256E-03 


80 


0.11989E+01 


0.80502E+02 


0.79970E4 


-02 


0.62812E-02 


0.37011E-03 


85 


0.12756E+01 


0.85501E+02 


0.84970E4 


-02 


0.58974E-G2 


0.35594E-03 


90 


0.13584E+01 


0.90501E+02 


0.89969E4 


-02 


0.55623E-02 


0.34087E-03 


95 


0.14472E+01 


0.95500E+02 


0.94969E4 


-02 


0.52652E-02 


0.32589E-03 


100 


0.15416E+01 


0.10050E+03 


0.99969E+02 


0.49994E-02 


0.31150E-03 



Table 2; Illustration of Theorems [^ \10\ with c = 100. The quantities Abovefn) 
and Below(n) are defined by (j253p . 



n 




Above(ri) 


Below(n) 


Abovc(?T,)— n 


n— Bclow(n) 


n 


n 


637 


0.10005EH 


-01 


0.63759E4-03 


0.63697E4 


-03 


0.93059E-03 


0.51797E-04 


640 


0.10025EH 


-01 


0.64055E-I-03 


0.63997E-f03 


0.85557E-03 


0.49251E-04 


645 


0.10063EH 


-01 


0.64552E-I-03 


0.64497E4 


-03 


0.80101E-03 


0.39996E-04 


650 


0.10105EH 


-01 


0.65051E-H03 


0.64997E4 


-03 


0.78412E-03 


0.39578E-04 


655 


0.10149EH 


-01 


0.65551E-h03 


0.65497E4 


-03 


0.77352E-03 


0.40527E-04 


660 


0.10195EH 


-01 


0.66050E-f03 


0.65997E4 


-03 


0.76512E-03 


0.41359E-04 


665 


0.10243E+01 


0.66550E+03 


0.66497E4 


-03 


0.75777E-03 


0.41942E-04 


670 


0.10292EH 


-01 


0.67050E-h03 


0.66997E-H03 


0.75103E-03 


0.42321E-04 


675 


0.10343EH 


-01 


0.67550E-h03 


0.67497E4 


-03 


0.74469E-03 


0.42547E-04 



Table 3: Illustration of TheoremsWl I J 01 with c — 1000. The quantities Above(n) 
and Below (n) are defined by (|253p . 



n 


Xn/C^ 


Above(n) 


Below(7i) 


Abovc{n) — n 


n— Below(n) 


n 


n 


1 


0.29824E-01 





10395E+01 


O.lOOOOE-fOl 


0.39511E-01 


0.00000E4-00 


9 


0.18531E-f00 





90625E+01 


0.89818E-I-01 


0.69444E-02 


0.20214E-02 


19 


0.36985E-I-00 





19069E4-02 


0.18981E-I-02 


0.36421E-02 


0.10180E-02 


29 


0.54240E-K00 





29075E-t-02 


0.28980E+02 


0.25825E-02 


0.69027E-03 


39 


0.70125E-H00 





39082E-H02 


0.38979E-t-02 


0.21102E-02 


0.53327E-03 


49 


0.84356E-H00 





49096E-K02 


0.48978E-t-02 


0.19543E-02 


0.45122E-03 


54 


0.90685E-t-00 





54110E-t-02 


0.53977E-t-02 


0.20330E-02 


0.43263E-03 


59 


0.96278E-K00 





59146E-K02 


0.58974E-I-02 


0.24725E-02 


0.44189E-03 


63 


0.99867E-F00 





63420E-t-02 


0.62966E-f02 


0.66661E-02 


0.53355E-03 



Table 4: Illustration of Theorems\^ \1(A with c = 100. The quantities Abovefn) 
and Below(n) are defined by (|254[) . 
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In Tabic m we illustrate Theorems [31 [TUlin the case Xn < with c — 100. 
The structure of Table S] is the same as that of Tables [TJ [5J [3] with the only 
difference: the third and fourth column contain the upper and lower bound on 
n given, respectively, via (|107l) in Theorems [5] and (|157l) in Theorem [TUl i.e. 



Below(n) 



Above(n) = — 





2 








TT 


2 




TT 


Jo 



1 -i2 



dt ~ 1 -\ \/Xn ■ E ( asin {tr, 

TT 



2 / c 



(254) 



where t„ and T are the maximal roots of ^pn and in the interval (—1,1), 
respectively. The values in the first row grow up to [2c/7rJ, in correspondence 
with Theorem [15] in Section 14.1.21 We observe that both bounds in the third 
and fourth columns are correct and the lower bound is always more accurate. 
This behavior is similar to that observed in Tables [1] [21 [3l 



n 


(n - 2c/7r- 1) /c 


Xn 


(f (n+1))' 


(f (n + l))Vxn-l 


640 


0.23802E-02 


0.10025E- 


h07 


0.10138E+07 


0.11248E-01 


660 


0.22380E-01 


0.10195E- 


h07 


0.10781E+07 


0.57443E-01 


680 


0.42380E-01 


0.10395E^ 


h07 


0.11443E+07 


0.10082E+00 


700 


0.62380E-01 


0.10615E^ 


h07 


0.12125E+07 


0.14229E+00 


720 


0.82380E-01 


0.10850E^ 


h07 


0.12827E+07 


0.18215E+00 


740 


0.10238E+00 


O.lllOOE^ 


h07 


0.13548E+07 


0.22054E+00 


760 


0.12238E+00 


0.11363E- 


h07 


0.14289E+07 


0.25757E+00 


780 


0.14238E+00 


0.11637E+07 


0.15050E+07 


0.29330E+00 


800 


0.16238E+00 


0.11923E+07 


0.15831E+07 


0.32777E+00 


820 


0.18238E+00 


0.12219E- 


h07 


0.16631E+07 


0.36105E+00 



Table 5: Illustration of Theorem \lS\ with c = 1000. 



n 


(n- 2c/7r- 1) /c 


Xn 


(f (n+1))' 


(f (n+l))Vx«-l 


6400 


0.32802E-02 





10022EH 


h09 





lOllOEH 


h09 


0.87670E-02 


6600 


0.23280E-01 





10191EH 


h09 





10751EH 


h09 


0.55007E-01 


6800 


0.43280E-01 





10390EH 


h09 





11413EH 


h09 


0.98410E-01 


7000 


0.63280E-01 





10609EH 


h09 





12094EH 


h09 


0.13991E+00 


7200 


0.83280E-01 





10845EH 


h09 





12795EH 


h09 


0.17979E+00 


7400 


0.10328E+00 





11094EH 


h09 





13515E+09 


0.21821E+00 


7600 


0.12328E+00 





11357EH 


h09 





14255EH 


h09 


0.25526E+00 


7800 


0.14328E+00 





11631EH 


h09 





15016EH 


h09 


0.29102E+00 


8000 


0.16328E+00 





11916EH 


h09 





15795EH 


h09 


0.32552E+00 


8200 


0.18328E+00 





12213EH 


h09 





16595EH 


h09 


0.35883E+00 



Table 6: Illustration of Theorem \13\ with c — 10000. 



In Tables [3 [H we illustrate Theorem [H with c = 1000 and c = 10000, 
respectively. The first column contains the PSWF index n, which starts from 
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roughly 2c/7r and increases by steps of c/50. The second column displays the 
normalized distance d„ between n and {Icjix + 1), defined via the formula 



n — Icjix — 1 



(255) 



The third column contains x„. The fourth and fifth column contain the upper 
bound on defined in Theorem [T31 and the relative error of this bound, 
respectively. We observe that the bound is slightly better for c = 10000, if 
we keep dn fixed. On the other hand, for a fixed c, this bound deteriorates as 
n grows. In fact, starting from n k, {2/t: + 1/6) • c, this bound becomes even 
worse than ^ (this value is n = 825 for c = 1000 and n = 8254 for c = 10000). 
Since Theorem [T51 is a simplification of more accurate Theorems [TUl the latter 
observation is not surprising. Nevertheless, the high accuracy for n 2c/7r and 
the simplicity of the estimate make Theorem [T3l useful (see also Figure [3]). 



10^ 



10 



-1 



10" 



10" 



c=1000 



xxx>oo<xxxxxxx>oooooooooooooooooooooooooooooooooo<x>oo<x>ooooooooooooooo< 




X error(UI) 
> error(U2) 
□ error(U3) 



640 



660 



680 



700 



Figure 3: Relative error of upper hounds on Xn with c — 1000, on the logarithmic 
scale. The bounds are defined, respectively, via (j256p . (|257p . (|258p . 



In Figure [31 we illustrate Theorems |31 [SI [101 via comparing the relative 
accuracy of the corresponding upper bounds on Xn- More specifically, we choose 
c = 1000, and, for each integer 630 < n < 710, we evaluate numerically the 
following quantities. First, we compute Xn (see (IT5t in Section [^?T|) . Second, 
we compute the upper bound on Xm defined via the right-hand side of (|16p of 
Theorem [H in Section I^TTl namely, 

[/i(n) = +n • (71 + 1). (256) 
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Third, we compute the upper bound on Xn, defined via (jl73p of Theorem 1131 
namely, 



U, 



(257) 



Finally, we compute the upper bound on x„, defined via (|217p of Theorem [501 
namely. 



[/3(n)=c^. l + if — + — 



2c 



2c 



(258) 



where H is defined via (|216l) in Theorem EDI In Figure [31 we plot the relative 
errors of Ui{n), U2{n), Uj,{n) as functions of n, on the logarithmic scale. 

We observe that J7i(n) significantly overestimates Xn, and the relative ac- 
curacy of Ui{n) remains roughly the same for all 630 < n < 710. On the other 
hand, the relative accuracy of U2{n) is higher than that of Ui{n); however, it 
deteriorates as n grows: from below 0.01 for n < 640 to above 0.1 for n > 680 
(see also Table [51 above) . Finally, C/3 (n) displays much higher relative accuracy 
than both Ui{n) and U2{n): the relative accuracy of C/3(n) remains below 0.004 
for aU 630 < n < 710. 



10^ 



10" 



10 



10" 



-2 



10 



-4 



c=1000 



»000«0<><>00<»<XXX><X)0<X><X>000<><X>0000<X>^^ 



> 
> 

> 




X error(LI) 
> error(L2) 
° error(L3) 



640 



660 



680 



700 



Figure 4: Relative error of lower bounds on Xn with c — 1000, on the logarithmic 
scale. The bounds are defined, respectively, via (j259p . (|260p . (j26ip . 



In Figure [H we illustrate Theorems [H [151 IlHl via comparing the relative 
accuracy of the corresponding lower bounds on More specifically, we choose 
c = 1000, and, for each integer 630 < n < 710, we evaluate numerically the 
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following quantities. First, we compute Xn (see ([T5|) in Section [^?T]) . Second, 
we compute the lower bound on defined via the left-hand side of (1161) of 
Theorem d] in Section [^?T1 namely, 

Li(n) = n- (n-t- 1). (259) 

Third, we compute the trivial lower bound on Xn, established in Theorem 1151 
namely, 

L2{n) = c^. (260) 

Finally, we compute the lower bound on x„ , defined via (|217p of Theorem 
namely, 

L3(n) = c2-(l + i/(^-l)), (261) 

where H is defined via (j216p in Theorem [20l In Figure |4j we plot the relative 
errors of Li{n), L2(n)^ L'i{n) as functions of n, on the logarithmic scale. 

We observe that Li(n) significantly underestimates x„, and the relative ac- 
curacy of Li{n) remains roughly the same for all 630 < n < 710. Even the 
trivial lower bound 1^2 (n) = displays a higher relative accuracy, which, obvi- 
ously, deteriorates as n grows. Finally, L^iiji) is much more accurate than both 
Li(n) and L2{n): the relative accuracy of L-i{n) remains below 0.0006 for all 
630 < n < 710. We also observe, that the relative accuracy of L-i{n) is about 
an order of magnitude higher than that of Uz{n), defined via (|258p above (see 
Figure El). 



i 


ti+l — ti 


/(t.+l)+i'(ti+l)/2 


TT 

fiu) 


lower error 


upper error 


44 
46 


0.27468E-01 
0.27453E-01 


0.27464E-01 
0.27439E-01 


0.27470E-01 
0.27460E-01 


0.13152E-03 
0.52432E-03 


0.63357E-04 
0.24265E-03 


60 
62 


0.26685E-01 
0.26437E-01 


0.26573E-01 
0.26303E-01 


0.26741E-01 
0.26506E-01 


0.42160E-02 
0.50867E-02 


0.21008E-02 
0.25968E-02 


70 
72 


0.24700E-01 
0.23948E-01 


0.24418E-01 
0.23602E-01 


0.24863E-01 
0.24158E-01 


0.11404E-01 
0.14473E-01 


0.66360E-02 
0.87772E-02 


84 
86 


0.96757E-02 
0.39568E-02 


0.81279E-02 
0.22125E-02 


0.10948E-01 
0.55074E-02 


0.15996E+00 
0.44083E+00 


0.13147E-I-00 
0.39188E+00 



Table 7: Illustration of Theorem \ll\ with c — 100 and n — 87. 

In Tables H [l we illustrate Theorems [IH [H with c = 100, n = 87 and 
c = 1000, n = 670, respectively. The first column contains the index i of the 
ith root ti of 4'n inside (—1,1). The second column contains the difference be- 
tween two consecutive roots t^+i and ti. The third and fourth columns contain, 
respectively, the lower and upper bounds on this difference, given via (|162p in 
Theorem [TTl The last two columns contain the relative errors of these bounds. 
We observe that both estimates are fairly accurate when ti is far from 1, and 
the accuracy increases with c. The best relative accuracy is about 0.01% for 
c = 100 and 0.0001% for c = 1000. Both bounds deteriorate as i grows to n. For 
both values of c the relative accuracy of the lower bound for i = n — 1 is as low 
as 44%, and that of the upper bound is about 39%. In general, the upper bound 
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i 


ti+i — U 


/(t>+i)+f(t, + i)/2 




lower error 


upper error 


336 
338 


0.30967E-02 
0.30967E-02 


0.30967E-02 
0.30967E-02 


0.30967E-02 
0.30967E-02 


0.19367E-05 
0.52185E-05 


0.59233E-06 
0.86461E-06 


400 
402 


0.30948E-02 
0.30947E-02 


0.30945E-02 
0.30944E-02 


0.30949E-02 
0.30947E-02 


0.11172E-03 
0.11547E-03 


0.10078E-04 
0.10427E-04 


500 
502 


0.30813E-02 
0.30808E-02 


0.30802E-02 
0.30797E-02 


0.30815E-02 
0.30810E-02 


0.37302E-03 
0.38101E-03 


0.41125E-04 
0.42311E-04 


601 
603 


0.30109E-02 
0.30071E-02 


0.30065E-02 
0.30025E-02 


0.30118E-02 
0.30080E-02 


0.14549E-02 
0.15168E-02 


0.30734E-03 
0.32775E-03 


667 
669 


0.10176E-02 
0.41703E-03 


0.85504E-03 
0.23323E-03 


0.11505E-02 
0.58020E-03 


0.15973E+00 
0.44073E+00 


0.13065E+00 
0.39128E+00 



Table 8: Illustration of Theorem \ll\ with c — 1000 and n = 670. 



c = 200 

-2 1 ' ' ' 




n 



Figure 5: Illustration of Theorem \2S\ with c = 200. Here i„ is the maximal root 
of Tpn in ("Ij 1); while the lower and upper hounds are defined via (j262p . (|263p . 
respectively. 
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is always more accurate. We also note that t^+i — U decreases monotonically as 
i grows, which confirms Theorem 1161 since Xn > (? in both cases. 

We illustrate Theorem E51 in Figure [5] We choose c — 200, and, for each 
integer 130 < n < 230, we evaluate numerically the following quantities. First, 
we compute the maximal root <„ of ipn in (—1,1). Second, we evaluate the 
eigenvalue Xn (see in Section [^TT|) . Then, we compute the lower and upper 
bounds on 1 — t„ , established in Theorem I23[ namely, 

TT^ 1 

lower (n) = ^ -, (262) 

^ ^ 8-(i + V2) xn-c2' y ) 

27r2 

upper{n) — -. (263) 

Xn - 

In Figure [5l we plot \og{lower(n)), \og{upper{n)) and log(l — i„), as functions 
of n. 

We observe that neither of (|262p . (|263p is a very accurate estimate of 1 — i„. 
Nevertheless, they correctly capture the behavior of 1 — 1„, up to a multiplicative 
constant. In particular, for all integer 130 < n < 230, 

l-t„ = ^M^, (264) 



where £,{n) is a real number in the range 



^2 



■{1 + V2) 



< ^{n) < 27^^ (265) 



as expected from Theorem [221 In other words, 1 — t„ is proportional to {xn 
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